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ABSTRACT 

We present new BeppoSAX LECS, MECS, and PDS observations of four flat-spectrum radio quasars 
(FSRQ) having effective spectral indices a ro and a ox typical of high-energy peaked BL Lacs. Our 
sources have X-ray-to-radio flux ratios on average ~ 70 times larger than "classical" FSRQ and lie 
at the extreme end of the FSRQ X-ray-to-radio flux ratio distribution. The collected data cover the 
energy range 0.1 — 10 keV (observer's frame), reaching ~ 100 keV for one object. The BeppoSAX band 
in one of our sources, RGB J1629+4008, is dominated by synchrotron emission peaking at ~ 2 x 10 16 
Hz, as also shown by its steep (energy index a x ~ 1.5) spectrum. This makes this object the first known 
FSRQ whose X-ray emission is not due to inverse Compton radiation. Two other sources display a 
flat BeppoSAX spectrum (a x ~ 0.7), with weak indications of steepening at low X-ray energies. The 
combination of BeppoSAX and ROSAT observations, (non-simultaneous) multifrcquency data, and a 
synchrotron inverse Compton model suggest synchrotron peak frequencies ss 10 15 Hz, although a better 
coverage of their spectral energy distributions is needed to provide firmer values. If confirmed, these 
values would be typical of "intermediate" BL Lacs for which the synchrotron and inverse Compton 
components overlap in the BeppoSAX band. Our sources, although firmly in the radio-loud regime, 
have powers more typical of high-energy peaked BL Lacs than of FSRQ, and indeed their radio powers 
put them near the low-luminosity end of the FSRQ luminosity function. We discuss this in terms of an 
anti-correlation between synchrotron peak frequency and total power, based on physical arguments, and 
also as possibly due to a selection effect. 

Subject headings: radiation mechanisms: non-thermal — galaxies: active, quasars — X-rays: galaxies 



1. INTRODUCTION 

Blazars constitute one of the most extreme classes of 
active galactic nuclei (AGN), distinguished by their high 
luminosity, rapid variability, high (> 3%) optical polar- 
ization, radio core-dominance, and apparent superlumi- 
nal speeds (Kollgaard 1994; Urry & Padovani 1995). The 
broad-band emission in these objects, which extends from 
the radio to the gamma-ray band, appears to be dom- 
inated by non-thermal processes from the heart of the 
AGN, often undiluted by the thermal emission present in 
other AGN. Therefore, blazars represent the ideal class to 
study to further our understanding of non-thermal emis- 
sion in AGN. 

The blazar class includes BL Lacertae objects, charac- 
terized by an almost complete lack of emission lines, and 
a subclass of radio quasars (which by definition display 
broad emission lines) which have been variously called 
highly polarized quasars (HPQ), optically violently vari- 
able quasars (OVV), core-dominated quasars (CDQ). One 



of their observational properties which is easier to define 
is the flat-spectrum radio emission {pt T < 0.5) and so we 
will refer to them as flat-spectrum radio quasars (FSRQ). 

Given the lack of prominent emission lines in BL Lacs, 
more than 95% of all known such objects have been discov- 
ered either in radio or X-ray surveys. Follow-up work on 
radio- and X-ray-selected samples has shown that the two 
selection methods yield objects with somewhat different 
properties. The energy output of most radio selected BL 
Lacs peaks in the IR/optical band (Giommi & Padovani 
1994; Padovani & Giommi 1995, 1996); such objects are 
now referred to as LBL (low-energy peaked BL Lacs). By 
contrast, the energy output of most X-ray selected BL 
Lacs (referred to as HBL: high-energy peaked BL Lacs) 
peaks at UV/X-ray energies. 

Padovani & Giommi (1995) and Sambruna, Maraschi 
& Urry (1996) have demonstrated that the difference in 
broad-band peaks for HBL and LBL is not simply phe- 
nomenological. Rather, it represents a fundamental dif- 
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ference between the two sub-classes. The location of the 
broadband peaks also suggests a different origin for the X- 
ray emission of the two classes. Namely, an extension of 
the synchrotron emission likely responsible for the lower 
energy continuum in HBL, which typically display steep 
(energy index a K <~ 1.5) X-ray spectra, and inverse Comp- 
ton emission in LBL, which have harder (a x ~ 1) spectra 
(Perlman et al. 1996; Urry et al. 1996; Padovani & Giommi 
1996). BeppoSAX observations of BL Lacs are confirming 
this picture (Wolter et al. 1998; Padovani et al. 2001; Beck- 
mann et al. 2002). 

In this respect, one could expect to find a similar range 
in peak frequencies in the FSRQ class - for which, until 
recently, no evidence existed. Indeed, it was suggested 
by some authors (Sambruna, Maraschi & Urry 1996), 
based upon the similarities of the optical-X-ray broad- 
band spectral characteristics of LBL and FSRQ, that no 
FSRQ with synchrotron peak emission in the UV/X-ray 
band should exist. 

Two studies have drastically changed this picture: 1. 
the multifrequency catalog of Padovani, Giommi, & Fiore 
(1997) identified 'more than 50 FSRQ (~ 17% of the 
FSRQ in their catalog) spilling into the region of param- 
eter space once exclusively populated by HBL; 2. about 
30% of FSRQ found in the deep X-ray radio blazar survey 
(DXRBS) (Perlman et al. 1998; Landt et al. 2001) were 
found to have X-ray-to-radio luminosity ratios, £ x /£ r , 
typical of HBL (L x /L r > 10~ 6 or a rx < 0.78), but 
broad (FWHM > 2,000 ^8-*) and luminous (L > 
10 43 erg s -1 ) emission lines typical of FSRQ. The discov- 
ery of a large population of "X-ray strong" FSRQ (la- 
beled HFSRQ by Perlman et al. (1998) to parallel the HBL 
moniker) represents a fundamental change in our percep- 
tion of the broadband emission of FSRQ. 

X-ray observations of these objects play a fundamen- 
tal role in finding their place within the blazar class. For 
example, if the X-ray spectra were found to be relatively 
steep, one could infer a dominance of synchrotron emis- 
sion, as observed in HBL. Flatter X-ray spectra, with 
corroborating evidence from the whole broad-band emis- 
sion (Padovani & Giommi 1996), would instead suggest 
inverse Compton emission. In the latter case, the simple 
equations LFSRQ = LBL (where by LFSRQ we mean the 
"typical" FSRQ with low-energy synchrotron peak) and 
HFSRQ = HBL would not be valid, and some more com- 
plicated explanation for the existence of this class should 
be sought. 

Our previous knowledge of the X-ray spectra of this new 
class of objects is scanty and mostly based on low signal- 
to-noise ratio (S/N) ROSAT data (Padovani et al. 1997; 
Padovani et al., in preparation) and is therefore also lim- 
ited to the relatively narrow 0.1 — 2.4 keV band. The Bep- 
poSAX satellite (Boella et al. 1997a), with its broad-band 
X-ray (0.1 — 300 keV) spectral capabilities, is particularly 
well suited for a detailed analysis of the individual X-ray 
spectra of these sources. 

In this paper we present BeppoSAX observations of four 
HFSRQ candidates, selected as described below. In § 2 
we present our sample, § 3 discusses the observations and 
the data analysis, while § 4 describes the results of our 
spectral fits to the BeppoSAX data. § 5 discusses the 
ROSAT data, § 6 presents the spectral energy distribu- 



tions and synchrotron-inverse Compton fits to the data, 
§ 7 discusses our results, while § 8 summarizes our conclu- 
sions. Throughout this paper spectral indices are written 
S v oc v~ a and the values H = 50 km s _1 Mpc -1 and 
qa = have been adopted. 

Readers wanting to skip the details of the data reduction 
and go directly to our results can read § 2, the summaries 
of our BeppoSAX and ROSAT results'in §§ 4.5 and 5.4, 
and then go straight to § 6. 

2. THE SAMPLE 

The sample selection was done in two separate steps, as 
our objects were observed in the BeppoSAX Cycle 2 and 
3. For both cycles we defined as HFSRQ flat-spectrum 
quasars with a ro , a ox , and a rx values within 2er from the 
mean values of the HBL in the multifrequency AGN cat- 
alog of Padovani, Giommi, & Fiore (1997). These are the 
usual effective spectral indices defined between the rest- 
frame frequencies of 5 GHz, 5,000 A, and 1 keV. X-ray 
and optical fluxes have been corrected for Galactic absorp- 
tion. The effective spectral indices have been K-corrected 
using the appropriate radio and X-ray spectral indices, 
available for most sources (otherwise mean values were as- 
sumed). Optical fluxes at 5,000 A have been derived and 
K-corrected as described in Padovani et al. (in prepara- 
tion) . We have chosen this definition so that our selection 
of HFSRQ matches the traditional "HBL box" as detailed 
in, e.g., Padovani & Giommi (1995). 




Fig. 1. — The a TO — a ox plane for flat— spectrum radio quasars 
(FSRQ) and HBL BL Lacs. Triangles represent FSRQ, crosses rep- 
resent HBL, while filled circles represent our sources. The region in 
the plane within 2cr from the mean o ro , o x, and a rx values of HBL 
is indicated by the solid lines. Data from the multifrequency AGN 
catalogue of Padovani, Giommi, & Fiore (1997). 

In January 1998, for Cycle 2, we selected the four X- 
ray brightest HFSRQ candidates then known. Three of 
these came from the AGN catalog of Padovani, Giommi, 
& Fiore (1997), while the fourth one was the X-ray bright- 
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Table 1 
Sample Properties. 



Name 


RA(J2000) 


Dcc(J2000) 


z 


O 


E 


FsGHz 


a r 


f>ro 


Oox 


Q!rx 


Galactic Nh 














mjy 








10 20 cm- 2 




WGA J0546.6-6415 


05 46 41.8 


-64 15 22 


0.323 


16.0 a 


14.7 b 


287 


-0.67 


0.39 


1.17 


0.66 


4.54 


RGB J1629+4008 


16 29 01.3 


+40 08 00 


0.272 


17.6 


17.1 


20 


-0.29 


0.34 


0.93 


0.54 


0.85 


RGB J1722+2436 


17 22 41.2 


+24 36 19 


0.175 


16.5 


15.1 


35 


-0.75 


0.26 


1.34 


0.63 


4.95 


S5 2116+81 


21 14 01.2 


+82 04 48 


0.084 


14.6 


13.3 


376 


0.26 


0.31 


1.41 


0.69 


7.41 



a J magnitude from the US Naval Observatory (USNO) A2.0 catalogue 
b F magnitude from USNO A2.0 catalogue 



Table 2 

BeppoSAX Journal of observations. 



Name 


LEGS 


LEGS 


MEGS 


MEGS 


PDS 


PDS 


Observing Date 




cxp. (s) 


count ratc a (cts/s) 


cxp. (s) 


count ratc a (cts/s) 


cxp. (s) 


count ratc a (cts/s) 




WGA J0546.6-6415 


18564 


0.036 + 0.002 


47234 


0.053 + 0.001 


20004 


0.148 ± 0.059 


1998 Oct 1-2 


RGB J1629+4008 


20512 


0.033 ± 0.002 


44759 


0.023 ± 0.001 


22798 


0.086 ± 0.054 


1999 Aug 11-12 


RGB J1722+2436 


12175 


0.008 ± 0.003 


43993 


0.053 + 0.001 


19732 




2000 Feb 13-14 


S5 2116+81 


13171 


0.147 + 0.004 


28826 


0.218 + 0.003 


13296 


0.286 + 0.073 


1998 Apr 29 


S5 2116+81 


5575 


0.107 + 0.006 


19674 


0.162 + 0.003 


10488 


0.240 + 0.084 


1998 Oct 12-13 



"Net count rate full band 



est HFSRQ in the DXRBS list of Pcrlman et al. (1998). 
We were granted BeppoSAX observing time for two of 
these sources. In January 1999, for Cycle 3, we were 
able to select even more extreme candidates by using 
the ROSAT-Green Bank survey (RGB) (Brinkmann et al. 
1997; Laurent-Muchlciscn ct al. 1998). Being based on 
the ROSAT All-Sky Survey (RASS) the RGB sample has 
an X-ray flux limit higher than DXRBS (/ x ~ 2 x KT 13 
erg cm -2 s -1 vs. ~ a few x 10~ 14 erg cm~ 2 s" 1 ) while it 
reaches slightly lower 5 GHz radio fluxes (/ r ~ 25 mjy vs. 
~ 50 mjy). RGB sources have then, by selection, higher 
fx/ fr ratios, that is lower a rx , as needed to get more ex- 
treme (with synchrotron peaks, ^ poa k, at higher energies) 
HFSRQ. The identification of FSRQ from the RGB sample 
is based on a cross-correlation with the NRAO VLA (Very 
Large Array) Sky Survey (NVSS) (Condon et al. 1998) 
at 1.4 GHz and is described in detail in Padovani et al. 
(in preparation). Amongst all FSRQ in the HBL region 
with a rx < 0.65 which, based on HBL spectral energy 
distributions (SEDs), would correspond to z/ pca k > 10 16 
Hz (Fossati et al. 1998), we again selected the four X- 
ray brightest HFSRQ candidates and, again, were granted 
BeppoSAX observing time for two out of four sources. 

The object list and basic characteristics are given in Ta- 
ble 1, which presents the source name, position, redshift, O 
(blue) and E (red) magnitudes from the Automatic Plate 
Measuring (APM) (Irwin et al. 1994), 5 GHz radio flux, 
radio spectral index, a ro , a ox , and a rx values, and Galactic 
N H ._ 

Figure 1 shows the a ro — a ox plane for FSRQ and HBL 
and the location of our sources, well within the region 
which includes ~ 95% of HBL and clearly offset from 
where most FSRQ are located. Note that, in fact, our ob- 
jects have (a rx ) = 0.63 ± 0.03. For comparison, the FSRQ 
with X-ray data belonging to the 2 Jy sample (~ 80% 
of the total), which include all "classical" FSRQ (e.g., 3C 
273, 3C 279, 3C 345, 3C 454.3: Padovani (1992)), have 
(ttrx) = 0.87 ± 0.01, which corresponds to an X-ray-to- 
radio flux ratio ~ 70 times smaller. We anticipate that 



our sources have powers which, although within the FSRQ 
range, are relatively low, as discussed in §§ 7.2.4 and 7.3. 
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Fig. 2.— The light curve of the two merged MECS units for RGB 
J1629+4008, along with its hardness ratio, defined as the ratio be- 
tween the count rates in the 4-10 kcV and 1.5-4 keV bands. 

3. OBSERVATIONS AND DATA ANALYSIS 

A complete description of the BeppoSAX mission is 
given by Boclla et al. (1997a). The relevant instruments 
for our observations arc the coaligncd Narrow Field Instru- 
ments (NFI), which include one Low Energy Concentra- 
tor Spectrometer (LECS) (Parmar et al. 1997) sensitive in 
the 0.1 - 10 keV band; three (two after May 1997) identi- 
cal Medium Energy Concentrator Spectrometers (MECS) 
(Boella et al. 1997b), covering the 1.5 - 10 keV band; 
and the Phoswich Detector System (PDS) (Frontera et al. 
1997), sensitive in the 15 — 300 keV band, coaligned with 
the LECS and the MECS. A journal of the observations is 
given in Table 2. 
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Table 3 

Single power-law fits, LECS + MECS. 



Name 






FlkeV 


^[2-10] 


^[0.1-2.41 


Norm 


xt/dof 


F-tcst, notes 




10 20 cm- 2 






erg cm -2 s' 1 


erg cm s 


(L/M) 




fixed-free Nh 


WGA J0546. 6-6415 


4.54 fixed 


0.72 ± 0.08 


0.64 ± 0.07 


3.84c-12 


4.14c-12 


0.82 


1.07/41 








0.72 ± 0.10 


0.64 ±0.08 


3.84c-12 


4.15e-12 


0.82 


1.09/40 


4% 


KOrS Jloz9+4008 


not; c „ ,j 

U.oo nxea 


1.5U ± 0.06 


n _l_ one; 
0.66 ± 0.05 


1.25e-12 


/ .96C-12 


0.67 


n on /oo 
0.90/28 






0.58l° :» 


1.45 ± 0.11 


0.63 ±0.06 


1.27c-12 


7.19c-12 


0.67 


0.88/27 


79% 


RGB J1722+2436 


4.95 fixed 


0.62 ± 0.22 


0.16 ±0.05 


1.07c-12 


9.77c-13 


0.81 


0.63/18 






1.02 (< 14.6) 


0.58l°; 2 3 


0.14 ±0.05 


1.08c-12 


8.85c-13 


0.79 


0.55/17 


92% 


S5 2116+81 


7.41 fixed 


0.73 ± 0.04 


2.71 ±0.17 


1.61e-ll 


1.77e-ll 


0.77 


0.96/61 


1998 Apr 29 




io.3t 3 ; 6 


0.75 ± 0.07 


2.83 ±0.22 


1.60e-ll 


1.87c-ll 


0.78 


0.91/60 


96% 




7.41 fixed 


0.77 ±0.07 


2.15 ±0.19 


1.19c-ll 


1.43e-ll 


0.67 


0.98/46 


1998 Oct 12-13 




7 oc+7.20 


0.77 ±0.09 


2 16+ - 22 


1.19c-ll 


1.44e-ll 


0.67 


1.00/45 


16% 




7.41 fixed 


0.74 ±0.04 


2.77 ±0.16 






0.70 


1.01/63 


sum 




g 6 + 3 - 7 
».D_ 2 g 


0.76 ± 0.06 


2.87 ±0.20 






0.70 


0.97/62 


92% 



Note. — The errors are at 90% confidence level for one (with fixed Nh) and two parameters of interest. The fit for S5 2116±81 includes 
PDS data. 



Table 4 

Broken power-law fits, LECS + MECS. 

Name N H as a H E break F lkeV -F[2-io] Norm X^Jdoi F-test 

10 2 " cm' 2 kcY ^Jy erg cm~ 2 s~ 1 (L/M) 

WGA J0546.6-6415 8.70 (< 20) l.l+^ g 0.6+",; 2 2.4+5'J 3 , 0.9+g 2 3.91c-12 075 0.99/38 88%f 94% 

RGB J1629±4008 b 0.85 fixed 1.5(unc.) 1.7(unc.) 3.8(unc.) 0.7(unc.) 1.22e-12 0.67 0.94/26 32% 

RGB J1722+2436 b 4.95 fixed 1.88(> 0.02) 0.6t o j 0.9(< 4.7) 0.1 ± 0.1 1.08c-12 0.80 0.59/16 75% 

S5 2116±81 7.41 fixed 0.4±°;g 0.76±°'og l.l±oj 2.7 ±0.3 1.60c-ll 0.77 0.92/59 92%, 4/29/98 
"Compared to single power— law with free Nh- 

'Errors arc at 90% confidence level for one parameters of interest. 

Note. — Unless otherwise indicated, errors are at 90% confidence level for two parameters of interest. The values of F-test refer to the 
comparison with single power-laws plus fixed Galactic column densities. 



The data analysis was based on the linearized, cleaned 
event files obtained from the BeppoSAX Science Data Cen- 
ter (SDC) on-line archive (Giommi & Fiore 1997). The 
data from the two MECS instruments were merged in one 
single event file by the SDC, based on sky coordinates. 
The event file was then screened with a time filter given 
by SDC to exclude those intervals related to events with- 
out attitude solution (i.e., conversion from detector to sky 
coordinates; see Fiore et al. (1999)). This was done to 
avoid an artificial decrease in the flux. As recommended 
by the SDC, the channels 1 — 10 and above 4 keV for 
the LECS and 0-36 and 220 - 256 for the MECS were 
excluded from the spectral analysis, due to residual cali- 
bration uncertainties. 

The spectral analysis was performed using the matri- 
ces and blank-sky background files released in November 
1998 by the SDC, with the blank-sky files extracted in 
the same coordinate frame as the source file, as described 
in Padovani et al. (2001). Because of the importance of 
the band below 1 keV to assess the presence of extra ab- 
sorption or soft excess (indicative of a double power-law 
spectrum), we have also checked the LECS data for differ- 
ences in the cosmic background between local and blank- 
sky field observations, comparing spectra extracted from 
the same areas on the detector (namely, two circular re- 
gions outside the 10' radius central region, located at the 
opposite corners with respect to the two on-board radioac- 
tive calibration sources). 

We looked for time variability in every observation, bin- 



ning the data in intervals from 500 to 4000 s, with null 
results except for RGB J1629±4008. This source, in fact, 
clearly varied during the observation (x 2 -test probability 
of constancy < 10~ 7 , see Fig. 2), doubling its MECS flux 
in about ~ 15 hours, with a more rapid decay of a factor 
of ~ 4 in about 7 hours at the end of the observation. 
The same pattern is present in the single MECS units 
light curves, and the background showed no significant 
flux variations. In spite of the flux variation, no appre- 
ciable spectral changes seem to be present. The hardness 
ratio, defined as the ratio between the count rates in the 
4-10 keV and 1.5-4 keV bands, is constant (with a x 2 ~ 
test probability ~ 35%). We have also looked for possible 
differences extracting the spectra in the "high" and "low" 
states (defined as > or < 0.03 cts/s, respectively), finding 
no significant differences between the spectral indices. The 
LECS data also show evidence of variability, although less 
significant than the MECS one (x 2 -test probability of con- 
stancy ~ 5 — 10%, depending on the bin size). We stress 
that our only variable source is also the only one whose 
X-ray band is clearly dominated by synchrotron emission 
(see next section). 

4. SPECTRAL ANALYSIS 

The spectral analysis was performed with the XSPEC 
10.0 package. Using the program GRPPHA in HEAsoft, 
the spectra were rebinned with more than 20 counts in 
every new bin, using the rebinning files provided by the 
SDC. Various checks using different rebinning strategies 
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WGA J0546.6-6415 LECS+MECS 
Fit: single powerlaw + galactic abs. 




channel energy (keV) 



WGA J0546.6-6415 LECS+MECS+PDS 
Fit: single powerlaw + galactic abs. 




+■ 



channel energy (keV) 



Fig. 3. — Best fit to the WGA J0546.6— 6415 data: a single power-law with Galactic column density. On the right we also include the PDS 
data. The PDS flux is above the LECS+MECS extrapolation due to the likely contribution of several X-ray sources in the field of view. See 
text for details. 



have shown that our results are independent of the adopted 
rebinning within the uncertainties. The LECS/MECS nor- 
malization factor was left free to vary in the range 0.65-1.0, 
as suggested by the SDC (Fiore et al. 1999). 

We fitted the combined LECS and MECS data both 
with single and broken power-law models, with Galac- 
tic and free absorption. The absorbing column density 
was parameterized in terms of Nh, the HI column den- 
sity, with heavier elements fixed at solar abundances and 
cross sections taken from Morrison & McCammon (1983). 
The Galactic value was derived from the nh program at 
HEASARC (based on Dickey & Lockman (1990)). The 
Nh parameter was set free to vary for all sources to check 
for internal absorption and/or indications of a "soft ex- 
cess" . The main results of our single power-law fits are 
reported in Tab. 3. The F-test probability quantifies how 
significant is the decrease in v 2 due to the addition of 
a new parameter (free Nh). The errors quoted on the 
fit parameters are the 90% uncertainties for one and two 
interesting parameters, for Galactic and free Nh respec- 
tively. The errors on the 1 keV flux reflect the statistical 
errors only and not the model uncertainties. The results 
for the broken power-law fit to the data are presented in 
Tab. 4. as and an are the soft and hard energy indices 
respectively. In this case the F-test probability quantifies 
the decrease in x 2 due to the addition of two parameters 
(from a single power-law fit to a broken power-law fit). 
The errors quoted in this case are the 90% uncertainties 
for two interesting parameters. 

We now discuss the details and results of the analysis 
for each objects. 

4.1. WGA J0546.6- 6415 

This source is located near the Large Magellanic Cloud 
(LMC) at - 51' from LMC-X3. Towards the edges of 
the LECS and MECS images (in the region at right as- 
cension < 05 ^ 46 m ), apparently diffuse emission is visible, 



consistent with the emission of some sources in the LMC 
already catalogued through ROSAT PSPC observations 
(see Haberl & Pietsch (1999); hereafter HP99). The most 
intense source is at ~ 11' from the quasar, with a count 
rate ~ 20% that of the quasar in the MECS images. It 
is spatially coincident with the source [HP99]0049, identi- 
fied as a foreground star in HP99. Due to the presence of 
serendipitous sources we have used an extraction region of 
4' also for the LECS instrument. 

The data are well fitted by a single power-law model 
with Galactic absorption. The best fit plots and values 
are reported in Fig. 3 and Tab. 3. The obtained spectral 
index is flat (a x < 1), thus indicating an inverse Compton 
power-law origin of the emission, given the SED charac- 
teristics of this quasar (see § 6). This result at first glance 
does not seem to confirm the "HBL-like" nature of this ob- 
ject, suggested by its location in the a ro — a ox plane. We 
note, however, that a slightly better fit, even if not signifi- 
cantly improved, can be obtained with a broken power-law 
model, allowing for a column density higher than Galactic 
(N H - 8.7 x 10 20 cm 2 , F-test - 88 - 94%, where the two 
values refer to a comparison with a single power-law with 
free and Galactic Nh respectively; see Tab. 4). In this 
case, the steeper component below ~ 2 keV may indicate 
the presence of the tail of the synchrotron emission, as 
characteristic of "intermediate" BL Lacs (e.g., Giommi et 
al. 1999). Variability could also play a role in this (see 
§5.1). 

Significant flux was detected also in the PDS instrument 
up to ~ 60 keV, but with a value not compatible with that 
expected by a simple extrapolation of the LECS+MECS 
spectrum (see Fig. 3, right). Due to the high number 
of X-ray sources in the field of view, we consider con- 
tamination the most likely explanation for this detection. 
In fact, there are 9 [HP99] sources within 40' from the 
quasar, with PKS 0552—640, an object classified as AGN 
in NED, at 40', and LMC-X3, a well known and bright 
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Fig. 4.— Single power-law fits to the LECS+MECS data of RGB J1629+4008, with absorption fixed at the Galactic value. On the right 
we include also the PDS data. The PDS flux is above the LECS+MECS extrapolation due to the contribution of A2199. See text for details. 



X-ray binary source, at ~ 51'. Given the known X-ray 
spectrum of the latter (Haardt et al. 2001), we checked 
its possible contamination level using the published best 
fit model, after accounting for the off-axis efficiency of the 
instrument. Indeed, the hardest component, a power-law 
of index a x ~ 1.7 (Haardt et al. 2001), although steep, 
is quite strong (~ 10/xJy at 1 keV), and provides a non- 
negligible contribution up to ~ 35 — 40 keV. 

4.2. RGB J1629+4008 

This source is located 37' from the center of the well 
known cluster Abell (A) 2199. The X-ray emission from 
the hot gas of the cluster is clearly visible at the edges 
of the LECS and MECS images. To avoid contamination 
from the cluster, therefore, an extraction region of 6' has 
been used for the LECS. The best fit results are reported 
in Fig. 4 and Tab. 3. The spectrum is well fitted by 
a single power-law with Galactic values for the absorbing 
column density. The resulting spectral index is quite steep 
(a x ~ 1.5), as typical for HBL objects. This source there- 
fore confirms its "HBL-like" properties also in the X-ray 
spectrum. We stress that RGB J1629+4008 is the first 
FSRQ ever found with synchrotron emission extending all 
the way to the X-ray band. 

The PDS detected this source up to ~ 35 keV but the 
flux is more than an order of magnitude higher than ex- 
pected from the LECS+MECS extrapolation (see Fig. 4, 
right). This is likely due to the cluster emission. A2199isa 
very X-ray bright source (T « 5 keV, F [2 _io] = 1.2 x 10~ 10 
erg s _1 , emission integrated over a region of radius 20' ; 
Dc Grandi & Molendi (2002) and private communication) 
and, being at 37' from the quasar, is entirely within the 
PDS field of view and near where the PDS efficiency is 
still ~ 50% (42' ). Although its emission decreases ex- 
ponentially towards higher energies, its high brightness 
still gives a flux in the 15-40 keV band ~ 4.4 x 10~ 12 
erg s _1 , higher than the quasar LECS+MECS extrapola- 
tion (~ 3.2 x 10 -13 erg s _1 ). As a check, we therefore 
fitted the LECS+MECS+PDS data including in the PDS 
the best fit model for the A2199 spectrum (De Grandi & 



Molendi 2002), after accounting for the off-axis efficiency. 
In such a case the PDS data are well accounted for, with 
a total xt~ 1-03 (as compared to 1.38 without the cluster 
contribution) . 

Note that for completeness we also fitted the data for 
this source with a broken power-law with absorption fixed 
at the Galactic value (see Tab. 4). No improvement with 
respect to the single power-law fit case was obtained. 

RGB J1722+243 LECS+MECS 
Fit: single powerlaw + galactic abs. 
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Fig. 5. — RGB J1722+243 single power-law model with Galactic 
absorption. 

4.3. RGB J 1722+2436 

This source is the weakest of the four observed in the 
soft X-ray range (see the LECS count rates in Tab. 2). 
To enhance the S/N, therefore, we used a 6' extraction re- 
gion for the LECS, as suggested by the SDC (Fiore et al. 
1999). The source is not detected by the PDS. The data 
are well fitted by a flat (a x < 1) single power-law model 
with Galactic absorption (see Fig. 5 and Tab. 3). The 
residuals, however, show the presence of some soft excess 
below 0.5—0.6 keV. This is confirmed by the free Nh fit, 
which gives values below the Galactic one. Due to the 
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Fig. 6. — BeppoSAX data and fitted spectra for our two observations of S5 2116+81. The data are fitted with a single power-law model 
with Galactic absorption. 



low S/N and the narrow band affected, the significance 
is only marginal (F-test ~ 92%, last column in Tab. 3). 
A fit with a broken power-law with the absorption fixed 
at the Galactic value yields a steep low-energy component 
with a x ~ 1.9, though unconstrained (see Tab. 4). The 
improvement with respect to the single power-law fit and 
Galactic absorption is however not significant (~ 75%) 

The flatness of the X-ray spectrum, as for WGA 
J0546. 6+6415, does not confirm the "HBL-like" nature 
of this objects as suggested by the broad band indices. 
However, the possible steeper component at soft X-ray 
energies, if confirmed, may be attributed to the tail of the 
synchrotron emission, thus locating the observed X-ray 
band in the "valley" between synchrotron and Compton 
emissions (see, in fact, Fig. 7). In this case, the quasar 
could be the FSRQ counterpart of "intermediate" BL Lacs. 
Again, variability could have also played a significant role. 

4.4. S5 2116+81 

This object was observed by BeppoSAX twice, six 
months apart (see Tab. 2). In both observations the source 
was detected in the PDS up to ~ 90 keV. The source varied 
between the two observations (~ 20%), with the highest 
state in the first observation. Both datasets are well fitted 
by a flat single power-law model (see Fig. 6 and Tab. 3), 
with very similar spectral indices (within the errors). The 
PDS data agree, within the errors, with the extrapolation 
of the LECS+MECS spectrum. Given the similarity of the 
spectral properties, we have also added the event files to 
obtain a higher S/N spectrum, whose fit is also reported 
in Tab. 3. 

The free Nh value agrees with the Galactic one in the 
second dataset (within the large errors), while there is a 
marginal indication of possible excess absorption in the 
first and the summed datasets, according to the F-test 
(96% and 92%, respectively). A broken power-law model 
with Galactic absorption provides an equivalent fit (giving 
roughly the same xt of the free Nh fit), with a flatter index 
(a x ~ 0.4) below 1.1 keV; however the significance of the 
improvement is slightly less (F-test ~ 92%; see Tab. 4), 



due to the higher number of parameters. The flat spectral 
indices in both observations (i.e., flux states) classify this 
source as a standard FSRQ. 

4.5. Summary 

To summarize our BeppoSAX spectral results, the fitted 
energy indices are fiat, a x < 0.8, for all but one source. 
The mean value is a x = 0.90 ± 0.20 and the weighted 
mean is a x = 1.03 ± 0.04. This latter value is heavily 
weighted towards RGB J1629+4008, the only object with 
a steep (a x ~ 1.5) slope. Excluding this source we de- 
rive a mean value a x = 0.70 ± 0.04 and a weighted mean 
a x = 0.75 ±0.05. 

Three of our sources may have a spectrum which is more 
complex than a single power-law, as also illustrated by the 
residuals to the single power-law fits. A double power-law 
model fitted to their spectra indicates, within the rather 
large errors, a low-energy (E < 1 keV) excess for WGA 
J0546.6+6415 and RGB J1722+2436, with a flatter com- 
ponent emerging at higher energies. In fact, the spec- 
tra are concave with quite a large spectral change, with 
c^s — a n ~ 0.5 — 1.3, and energy breaks around E ~ 1 — 2 
keV. The hard X-ray spectral index is an ~ 0.6, while 
as ~ 1.1 — 1.9. Although the fit is improved by using 
a double power-law model, an F-test shows that the im- 
provement is more suggestive than significant, with proba- 
bilities - 75-94%. We do include RGB J1722+2436 here, 
despite its relative low F-test probability, as we believe 
that a best fit N-r value a factor ~ 5 below the Galac- 
tic one is highly unlikely and strongly suggestive of a soft 
excess. We also note that by adding up the x 2 values in 
Tab. 4, an F-test shows that the improvement in the fit 
provided by a double power-law model for the two sources 
together is significant at the ~ 93% level. S5 2116+81, on 
the other hand, shows a hint of a flatter component at low 
energies (F < 1 keV), with as — oh ~ —0.4. 

In short, one of our sources appear to show an "HBL- 
like" X-ray spectrum, while the other three are character- 
ized by flat "LBL-like" spectra. One, or possibly two, of 
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Table 5 

WGA J0546. 6-6415. ROSAT Journal of observations. 



Obs. ID 


Instrument 


Off axis 


Exposure time 


Net count rate 


Observing date 






(arcmin) 


(sec.) 


(cts/s) a 




rpl40008n00 


PSPC-C 


37.63 


1334 


0.103 ±0.012 


1990 Jul 9 


rpl40638n00 


PSPC-C 


37.63 


1500 


0.101 ±0.011 


1990 Jul 10 


rpl40009n00 


PSPC-C 


21.98 


1543 


0.124 ±0.011 


1990 Jul 9 


rpl40639n00 


PSPC-C 


21.98 


1628 


0.140 ±0.011 


1990 Jul 10 


rp400078n00 


PSPC-B 


51.66 


7490 


0.156 ±0.005 


1993 Apr 10 



a Net count rate full band 



Table 6 

WGA J0546. 6-6415. ROSAT single power-law fits. 



Obs. 


N H 

10 2() cm- 2 




FlkeV 

^Jy 


^[0.1-2.4] 
erg cm~ 2 


Xt/dof 

s- 1 


PSPC-B 


4.54 fixed 


1.29 ± 0.10 


0.88 ±0.06 


8.69e-12 


0.61/31 


PSPC-C sum8 


4.54 fixed 




0.45 ± 0.05 


3.06c-12 


1.09/16 


PSPC-C sum9 


4.54 fixed 


o.s±8:I 


0.48 ± 0.04 


3.34e-12 


1.01/14 


PSPC-C fit total 


4.54 fixed 




0.47 ±0.03 


3.22e-12 


1.00/32 



Note. — The errors are at 90% confidence level for one parameter of interest. 



these sources, however, show hints of concave X-ray spec- 
tra typical of "intermediate" BL Lacs. 

5. ROSAT PSPC DATA 

In order to study possible variations of the X-ray spec- 
tral properties (and synchrotron peak) in these objects and 
have a better understanding of the underlying emission 
processes, we have also looked for other X-ray observa- 
tions available, in particular from the ROSAT archive, to 
take advantage of the higher resolution and collecting area 
at low energies. We found Position Sensitive Proportional 
Counter (PSPC) observations for two sources, namely 
WGA J0546.6-6415 and RGB J1629+4008. These ob- 
jects are located near two famous and well observed X- 
ray sources, i.e., LMC-X3 and A2199 respectively, and 
thanks to the wide field of view of the ROSAT satellite we 
could find observations where our targets were present as 
serendipitous sources. The journals of the ROSAT obser- 
vations are reported in Tab. 5 and Tab. 7. A third object, 
RGB J1722+2436, was detected in the RASS (Voges et al. 
1999). We now present the results of the ROSAT analysis, 
along with a comparison with the BeppoSAX results. 

5.1. WGA J0546.6- 6415 

Many PSPC observations of the Large Magellanic Cloud 
region are present in the ROSAT public archive, both for 
B and C detectors, but most have short exposures (< 1000 
s) and/or with the quasar at high off- axis angles. We 
therefore decided to analyze the best ones, i.e., the longest 
observation (7490 s, PSPC-B, pointed on LMC-X3), and 
those with the lowest off-axis displacement (PSPC-C). Ta- 
ble 5 reports the basic information, namely the observation 
sequence ID code, the detector used, the off-axis position 
of our quasar in the field of view, the exposure time, the 
net count rate full band and date of observation. The 
PSPC-C observations took place in the same epoch (1990 
July 9 and 10), with the satellite often changing position 
between the two pointings (thus the two off-axis positions 



in the same day). The many short observations have been 
then reduced and summed together in a single dataset by 
the ROSAT data center. 

We fitted all the PSPC-C datasets independently, with a 
single power-law model plus Galactic absorption, obtain- 
ing good fits with similar values for the spectral indices 
and normalizations, within the errors. To increase the 
S/N, therefore, we summed the event files for the two ob- 
servations with the same source detector coordinates (i.e., 
same off-axis angle). The fit results are reported in Tab. 
6, where "sum8" refers to the summed event files ..008n00 
and ..638n00, and "sum9" refers to the summed event files 
..009n00 and ..639n00. Given the similarity of the values, 
we have also fitted the two datasets together (sum8 and 
sum9), using the same model parameters for both (includ- 
ing normalization). The resulting fit is good (xt = 1-00) 
thus confirming that both datasets have compatible spec- 
tral properties. 

All datasets analyzed are well fitted by single power-law 
models with Galactic column density. The free Nh fits pro- 
vide values around the Galactic ones, and other absorbing 
models and/or a broken power-law model do not improve 
the fits. The results seem to indicate a variation of the 
spectral properties for this source: the 1993 PSPC-B spec- 
trum is characterized by a steep "HBL-like" slope (~ 1.3), 
contrary to what obtained by both the 1990 PSPC-C and 
1998 BeppoSAX observations, which agree with each other 
reporting a flat and roughly equal spectral index (~ 0.8 
and 0.7 respectively). The 1990 and 1998 flux levels are 
similar, too, with the BeppoSAX flux about 30% higher 
than the ROSAT one. Compared to these, instead, in 
1993 the quasar was in a high state, with a 0.1-2.4 keV 
flux more than 2.5 times higher than the 1990 values. The 
X-ray spectrum of this object, therefore, seems to have a 
"steeper when brighter" behavior, which may indicate a 
shift towards higher energies of the synchrotron peak dur- 
ing high states, when synchrotron emission dominates over 
inverse Compton in the X-ray band (Padovani & Giommi 
1996; Giommi et al. 1999). 
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Table 7 

RGB J1629+4008. ROSAT Journal of observations. 



Obs. ID 


Instrument 


off axis 
(arcmin) 


exposure time 
(sec.) 


net count rate 
(cts/s) a 


Observing date 


rpl50083n00 
rp800644n00 
rp701507n00 


PSPC-C 
PSPC-B 
PSPC-B 


35.28 
35.28 
00.29 


10553 
40999 
5175 


0.641 ±0.008 
0.676 ± 0.004 
0.960 ±0.014 


1990 Jul 18 
1993 Jul 25 
1993 Jul 30 



a Net count rate full band 



Table 8 

RGB J1629+4008. ROSAT single power-law fits. 



Obs. Nh ~* F lkeV F[ . 1-2.41 X^/dof 

10 20 cm -2 /jjy erg cm" s _1 



rpl50083n00 


0.85 fixed 


2.03 ±0.05 


0.33 ± 0.02 


8.15c-12 


2.16/53 




0.7 ±0.2 


1.95 ±0.14 


0.34 ± 0.03 


7.43c-12 


2.17/52 


rp800644n00 


0.85 fixed 


r, r,o + 0.02 

z - za -0.04 


0.29 ± 0.01 


9.68c-12 


2.12/52 




0.82 ±0.01 


9 91 +0.08 
z - zl -0.04 


0.30 ±0.01 


9.39e-12 


2.14/51 


rp701507n00 


0.85 fixed 


o ic+0.06 
z - ld -0.05 


0.33 ± 0.03 


9.48e-12 


1.37/43 




0.4 ±0.2 


i-88iu:i 4 4 


0.34 ±0.03 


6.74e-12 


0.95/42 



Note. — The errors are at 90% confidence level for one (fixed Nh) or two parameters of interest. 



It must be noted, however, that these considerations are 
based on the 1993 PSPC-B observation, which may be af- 
fected by possible calibration problems. In fact, it has 
the highest off-axis angle (it is near the edge of the field 
of view), and these regions are usually not well studied 
and calibrated. In addition, also possible PSPC-B spec- 
tral miscalibrations, as reported in Iwasawa, Fabian, & 
Nandra (1999) and Mineo et al. (2000), must be taken 
into account, as these might steepen the ROSAT slopes 
by ~ 0.2 - 0.3 (Mineo ct al. 2000). 

5.2. RGB J1629+4008 

As for the previous object, we looked in the ROSAT 
public archive for the best observations possible, i.e., the 
longest ones and/or those with the lower off- axis angle. 
We found two very long observations of the cluster A2199 
taken in 1990 and 1993 with the PSPC-C and B instru- 
ments, respectively. These included our target in the field 
of view with an offset of ~ 35'. We found also a pointed ob- 
servation, taken 5 days after the long exposure on A2199. 
The basic information for these datasets is reported in Tab. 
7. 

At first we fitted all datasets with a single power-law 
model, both with Galactic and free absorption, but these 
did not provide an adequate fit (see xl m Tab. 8). The 
presence of a more complex spectral shape is also in- 
dicated by the free Nh fit of the on-axis observation 
(rp701507n00), which provided a good xl onr y with ab- 
sorption values below the Galactic ones. 

We therefore tried alternative models, reported in Tab. 
9. The fits were significantly improved both adopting a 
broken power-law model with Galactic absorption (top 
three rows in Tab. 9), and adding a warm absorption 
component, at the redshift of the A2199 cluster, to the 
single power-law model (bottom three rows in Tab. 9). 
This under the hypothesis that, at ~ 35' from the cen- 
ter, there is still a not negligible amount of cluster gas, 



likely partially ionized. The model used for the warm ab- 
sorber is ZWNDABS in XSPEC, an approximation to the 
warm absorber using the photoelectric cross-sections in 
Balucinska-Church & McCammon (1992). 

The results are a bit puzzling and do not settle clearly 
the question of which model should be preferred. The F- 
test does not help here, given that the number of model 
parameters is the same. The off-axis observations seem to 
be better fitted by the warm absorber model, according 
to the x 2 -test (15-20% against 2%, see the xl in Tab. 9), 
but this is not confirmed by the on-axis observation, which 
has the best xl f° r the broken power-law model. More- 
over, the obtained values for the extra (warm) absorption 
are different (at the 90% level), a fact which would imply 
a variable absorption with timescales less than five days 
(the time interval between the last two observations, i.e., 
rp800644 and rp701507), quite unlikely. For this reason, 
and possible feature calibration problems for the off-axis 
observations (outside the PSPC central rib ring), we con- 
sider the broken power-law model a better representation 
of these data. 

In both cases, however, the spectral properties of the 
continuum for this object confirm its "HBL-like" char- 
acter showed in the BeppoSAX observations, displaying 
steep spectra in all three datasets. This result is con- 
firmed also accounting for the possible ROSAT miscali- 
bration (i.e., subtracting a maximum Aa ~ 0.2 — 0.3). 
The broken power-law fits, which show a flattening of the 
spectrum towards higher energies, may indicate that the 
inverse Compton contribution is becoming not negligible, 
i.e., that the region where the synchrotron and Compton 
emissions cross was only slightly beyond the ROSAT band. 

5.3. RGB J1722+2436 

This source is included in the RASS. Yuan et al. (1998) 
give a 0.1-2.4 keV flux ~ 4.2 x 10~ 12 erg cm -2 s _1 , i.e., a 
factor ~ 4 larger than the BeppoSAX value, with a spectral 
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Table 9 

RGB J1629+4008. Other fits to the ROSAT datasets. 



Name 


10 2() cm" 2 


as 


an 


keV 


FlkeV 


-^[0.1 — 2.4] 
erg cm~ 2 s — 1 


Xt/doi 


F-test 


rpl50083n00 


0.85 fixed 


9 o+U.!2 
z '°-0.1 


i o +u - 4 

'-"-0.5 


0.7 ±0.2 


0.23 ±0.05 


8.70e-12 


1.47/51 


> 99.9% 


rp800644n00 


0.85 fixed 


9 r,r, + 0.09 


1 4+0-4 
'-^-0.7 


0.8 ± 0.2 


0.25 ±0.02 


9.95c-12 


1.43/50 


> 99.9% 


rp701507n00 


0.85 fixed 


9 7+3.8 
z -'-0.4 


1 q+0-2 


3+ ' 2 


0.13 ±0.08 


1.08c-ll 


0.91/41 


> 99.9% 






N]-j warm 




^window 


FlheV 


^[0.1-2.41 


Xt/doi 


F-tcst 






z = 0.03 




keV 




erg cm - s -1 






rpl50083n00 


0.85 fixed 


11 ±4 


1.75 ± 0.09 




0.54 ±0.06 


8.77e-12 


1.17/51 


> 99.9% 


rp800644n00 


0.85 fixed 


6 5+ 3 ' 3 
D -°-2.4 


2.07 ±0.07 


49+0.07 
u -^ y -0.27 


0.38 ±0.03 


9.78c-12 


1.27/50 


> 99.9% 


rp701507n00 


0.85 fixed 


1 6+ 3 ' 5 
l.D_ 1 2 


9 9+0.2 

z - z -o.i 


95+ 03 


0.36 ±0.05 


1.10e-ll 


1.04/41 


>99.8% 



Note. — Unless otherwise indicated, errors are at 90% confidence level for two parameters of interest. The values of the F-test refer to the 
comparison with single power-laws plus fixed Galactic column densities. 



index derived from the hardness ratio a x = 1.27_ a28 , as- 
suming Galactic absorption. This is steeper than the Bep- 
poSAX value and again suggests a "steeper when brighter" 
behavior. 

5.4. Summary 

To summarize our ROSAT spectral results, the fitted 
energy indices are relatively steep, a x ~ 1.4, steeper than 
the BeppoSAX values by Aa « 0.5. Possible ROSAT mis- 
calibrations are supposed to account for a steepening of 
the order Aa <~ 0.2 — 0.3, so we regard this difference 
as significant. In two cases the source was also brighter, 
suggesting a "steeper when brighter" behavior, indicative 
of a shift of the synchrotron peak towards higher energies 
during high states. Overall, this indicates a non- negligible 
synchrotron component in the X-ray band, which actually 
dominates for RGB J1629+4008. 

6. SPECTRAL ENERGY DISTRIBUTIONS 

To address the relevance of our BeppoSAX data in terms 
of emission processes we have assembled multifrcquency 
data for all our sources. The main source of information 
was NED, and the data are not simultaneous with our 
BeppoSAX observations. We also looked for near-infrared 
observations in the Two Micron All Sky Survey (2MASS) 
Second Incremental Release (Cutri et al. 2000), finding 
data only for WGA J0546. 6-6415. 

6.1. Optical Variability 

Given that our fits are mostly based on X-ray and opti- 
cal data and that the latter are particularly important to 
determine ^peak, we need to address the non-simultaneity 
of the optical and X-ray observations. The optical fluxes 
of FSRQ are in fact known to show significant variability, 
which could influence our results. We have then investi- 
gated the range of optical variability of our sources, with 
the following results. 

• WGA J0546. 6-6415. This is a newly identified 
source from the DXRBS (Perlman et al. 1998). 
Nevertheless, we did find variability information 
using SuperCOSMOS Sky Survey (Hambly et al. 
2001) data, spanning the period 1975 - 1986. In 



particular, the UK Schmidt telescope (UKST) and 
European Southern Observatory (ESO) red obser- 
vations, taken at the end of 1985/beginning of 1986 
and separated by 50 days, agree within 0.01 mag- 
nitude, and give R ~ 15.2. This is 0.5 magni- 
tudes fainter than the F = 14.7 magnitude given 
in the USNO A2.0 catalogue, which refers to mid 
1980. The Guide Star Catalogue-II (GSC-II), avail- 
able at http://www-gsss.stsci.edu/, gives J = 
15.9 ± 0.2 and F = 15.3 ± 0.2 for 1989, which indi- 
cates a brightening of ~ 0.6 magnitudes compared 
to the UKST blue value of 16.5 for 1975. 

• RGB J1629+4008. This is a Faint Images of the 
Radio Sky at Twenty centimeters (FIRST) (Becker 
et al. 1995) quasar. Helfand et al. (2001) have 
studied the long-term optical variability of FIRST 
quasars. They found that in 1953 Palomar Observa- 
tory Sky Survey (POSS I) observations this object 
was ~ 0.4 magnitudes fainter in B and ~ 0.6 mag- 
nitudes fainter in R than in 1996 CCD observations. 
The GSC-II gives J = 18.0 ±0.4 and F = 17.5 ±0.4 
for 1993. This is perfectly consistent with the CCD 
B and R values given by Helfand et al. (2001). 

• RGB J1722+2436. This source was identified as 
a quasar by Bond et al. (1977) due to its optical 
variability. These authors quote a range of photo- 
graphic magnitude for this source between 15.7 and 
16.7, with long "still-stands" between brief periods 
of activity. The GSC-II gives J = 16.7 ± 0.4 and 
F = 15.2 ± 0.4 for 1996, in good agreement with 
the O and E APM values given in Tab. 1. 

• S5 2116+81. The GSC-II gives J = 14.4 ± 0.4 for 
1993, consistent with the APM O value given in 
Tab. 1. 

In summary, within the available optical data, our 
sources display maximum variability amplitude ~ 0.5 — 1 
magnitudes. We have then conservatively assigned to the 
optical fluxes an error of 0.2 dex, which corresponds to 
±0.5 magnitudes (see Fig. 7). 



BeppoSAX Observations of Synchrotron X-ray Emission from Radio Quasars 



11 



6.2. Synchrotron Inverse Compton Model Fits 

The SEDs for our sources are shown in Figure 7. The 
BeppoSAX data have been converted to vj v units using 
the XSPEC unfolded spectra after correcting for absorp- 
tion. ROSAT data are shown by a bow-tie that represents 
the spectral index range. 

To derive the intrinsic physical parameters which could 
account for the observed data we have fitted the SED of 
our sources with a homogeneous, one-zone synchrotron in- 
verse Compton model as developed in Ghiscllini, Celotti, 
& Costamante (2002). This model is very similar to the 
one described in detail in Spada et al. (2001; it is the 
"one-zone" version of it), and is characterized by a finite 
injection timescale, of the order of the light crossing time 
of the emitting region (as occurs, for example, in the in- 
ternal shocks scenario, where the dissipation takes place 
during the collision of two shells of fluid moving at dif- 
ferent speeds). In this model, the main emission comes 
from a single zone and a single population of electrons, 
with the particle energy distribution determined at the 
time ti n j, i.e., at the end of the injection, which is the time 
when the emitted luminosity is maximized. Details of the 
model can be found in Ghiscllini, Celotti, & Costamante 
(2002), who have applied it successfully to both low-power 
BL Lacs and powerful FSRQ. Here we summarize its main 
characteristics. 

The source is assumed cylindrical, of radius R and width 
AR 1 = R/T (in the comoving frame, where T is the bulk 
Lorentz factor). The particle distribution N(-~y) is assumed 
to have the slope n [^(7) oc 7""] above the random 
Lorentz factor j c , for which the radiative cooling time 
equals the injection time. The latter is assumed to be equal 
to the light crossing time (i.e., i; n j ~ AR/c). The electron 
distribution is assumed to cut-off abruptly at 7 max > 7c- 
Below 7 C there can be two cases, depending on the values 
of 7 C and 7 min : 

i) If 7c > 7min, we have N(j) oc 7"("^ 1 ) between 

7min 

and 7 C and AT (7) oc 7 _1 below 7 m i n . 

ii) Alternatively, if 7 C < 7min, then N(-f) oc 7 -2 between 
7 C and 7 m ; n and ^(7) oc 7" 1 below 7 C . 

According to these assumptions, the random Lorentz 
factor 7 pca k of the electrons emitting most of the radiation 
(i.e., emitting at the peaks of the SEDs) is determined by 
the relative importance of radiative losses and can assume 
values in the range 7 m ; n to 7 max . 

Our sources have broad lines and therefore the contri- 
bution of the disk to the SED could be relevant. Further- 
more, photons produced in the broad line region could 
contribute to the seed photon distribution for the in- 
verse Compton scattering. We accounted for this by 
assuming that a fraction <~ 10% of the disk luminosity 
Ldisk i s reprocessed into line emission by the broad line 
region (BLR), £blr, assumed to be located at .Rblr- 
Lblr was estimated following Celotti, Padovani, & Ghis- 
cllini (1997) from the fluxes of the strongest broad lines. 
These were available for all our sources from our own 
spectrum (WGA J0546.6-6415; Perlman et al. (1998)), 
from published spectra kindly made available to us by 
Sally Laurent-Muehleisen (RGB J1629+4008 and RGB 
J1722+2436) and from Stickel, Kiihr, & Fried (1993) (S5 
2116+81). i?BLR is assumed to scale as • , following 
Kaspi et al. (2000). Disk emission is assumed to be a 



simple black-body peaking at 10 15 Hz. Within these as- 
sumptions, we can fix two important parameters for our 
modeling, and wc can reliably estimate the importance of 
the seed photons produced externally to the jet for the 
Compton scattering process. 

The source is assumed to emit an intrinsic luminosity L' 
and to be observed with the viewing angle 9. The model 
parameters are listed in Table 10, which gives the name of 
the source in column (1), L' in column (2), Ldisk in column 
(3), -Rblr in column (4), the magnetic field B in column 
(5), the size of the region R in column (6), the Lorentz 
factor T in column (7), the viewing angle 9 in column (8), 
the slope of the particle distribution n in column (9), the 
minimum Lorentz factor of the injected electrons 7 m j n in 
column (10), and 7 poa k and f^k in columns (11) and (12) 
respectively. Column (13) gives the magnetic field B cq 
derived assuming equipartition of the energy content of 
the magnetic field with the electron energy. Note that L' , 
B, R, T, 9, 7 m ; n and n are varied during the fitting pro- 
cedure, idisk and Lblr are fixed by observations, while 
7 P eak, t'peak' an d are derived quantities. 

In the case of a pure synchrotron self-Compton model, 
all the above parameters are constrained in sources for 
which: 1) we have an estimate of the minimum timescale of 
variability; 2) both the synchrotron and the self-Compton 
peak are well defined; 3) the spectral slopes below and 
above the peaks are known; 4) the redshift is known. As 
discussed in Tavecchio, Maraschi, & Ghisellini (1998), this 
suffices to fix the values of the magnetic field, the intrin- 
sic power of the source, the slopes of the emitting elec- 
tron distribution, the relativistic Doppler factor, and the 
dimension of the source. When the radiation produced ex- 
ternally to the jet is important there is one unconstrained 
unknown, but the superluminal motion of the radio knots 
observed in blazars indicates values of the bulk Lorentz 
factor in the range 10-15 on average, and we therefore use 
these values for our fits (see, e.g., Ghiscllini et al. 1998). 

We do not have information about the high-energy 
(Compton) peak for the sources in our sample so we lack 
the direct determination of 7 pC ak- But in the model we 
use here the finite time of injection of particles plays a 
crucial role, and we have an additional constraint with re- 
spect to the simplest synchrotron inverse Compton model. 
Namely, the peak of the synchrotron emission is either due 
to the electrons injected with 7 m ; n , or it is due to the elec- 
trons which were able to radiatively cool in the timescale 
tinj (i.e., electrons with Lorentz factor 7 C ). In the latter 
case (which is verified for all the sources studied in this 
paper but RGB J1629+4008), the value of 7 C depends on 
the total energy density (magnetic plus radiative) as seen 
in the comoving frame. Since the "external" radiation en- 
ergy density is known once we specify the bulk Lorentz 
factor (it is assumed to be produced by the BLR), 7 C de- 
pends on the magnetic field. The break in the particle 
distribution at 7 C corresponds to a break of the emitted 
synchrotron spectrum (and can often correspond to the 
peak of the synchrotron spectrum) and in this way we can 
reliably estimate the value of the magnetic field. 

We are aware of the fact that for three of our objects 
the shape, sampling, and lack of simultaneity of the SED 
do not allow us to firmly constrain the peak of the syn- 
chrotron emission (the exception being RGB J1629+4008), 
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Fig. 7. — Spectral energy distributions of our sources. Open points represent data from the literature (NED) while the BeppoSAX data 
are indicated by filled points. ROSAT data are shown by a bow-tie that represents the spectral index range. The dashed lines correspond to 
the one-zone homogeneous synchrotron and inverse Compton model calculated as explained in the text, with the parameters listed in Table 
10. Dotted lines represent the disk emission, while solid lines indicate the sum of disk plus non-thermal emission. The errors on the optical 
fluxes correspond to ±0.5 magnitudes (see text for details). 
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also because of the non-negligible contribution of the ther- 
mal disk component to the optical flux (especially for 
WGA J0546.6-6415 and S5 2116+81; see Fig. 7). How- 
ever, some reasonable arguments can be made. For WGA 
J0546.6-6415 and RGB J1722+2436 the presence of a 
steeper (and variable) component at soft X-ray energies 
(§§ 4.5 and 5.4) may be attributed to the tail of the syn- 
chrotron emission (see also Sect. 7.2.3), thus suggesting 
a synchrotron peak which cannot be at too low energies 
(> io 14 Hz). This is corroborated also by the IR and 
optical data for RGB J1722+2436, whose optical flux is 
not dominated by the disk emission. In the case of WGA 
J0546.6— 6415, which has a stronger disk component, the 
synchrotron peak could reach ~ 10 13 Hz only if one ne- 
glects the evidence for an upturn at low X-ray energies. S5 
2116+81 is the less constrained source, but it is also the 
case for which we have no indication of a steep, soft X-ray 
component. The radio to IR data suggest a peak above 
10 13 Hz and no set of parameters was found to be compat- 
ible with the IR and BeppoSAX data and with a z^ak ^ ess 
than ~ a few x 10 13 Hz. For all three objects, however, 
the synchrotron peak frequency has to be < 10 16 Hz, since 
otherwise the X-ray band would have a strong and obvi- 
ous steep synchrotron emission (as for RGB J1629+4008), 
contrary to what observed. In summary, we have modelled 
our objects mainly on the basis of the BeppoSAX (and 
partly ROSAT) data, trying to fit the other data as far as 
they could be compatible with the BeppoSAX ones. The 
derived parameters, then, apart from RGB J1629+4008, 
should be regarded as tentative. 

The model fits are shown in Fig. 7 as solid lines. The 
applied model is aimed at reproducing the spectrum orig- 
inating in a limited part of the jet, thought to be respon- 
sible for most of the emission. This region is necessar- 
ily compact, since it must account for the fast variability 
shown by all blazars, especially at high frequencies. The 
radio emission from this compact region is strongly self- 
absorbed, and the model cannot account for the observed 
radio flux. This explains why the radio data are system- 
atically above the model fits in the figure. Our model fails 
also to reproduce the optical flux of S5 2116+81. This 
is the source with the lowest redshift (z = 0.084) and we 
therefore expect a non-negligible contamination from the 
host galaxy, which is clearly detected in the Digitized Sky 
Survey (DSS) plates. 

Note that the model fits predict that the high energy 
inverse Compton emission has a luminosity comparable to 
the synchrotron one in all sources but one, S5 2116+81. 
The shape of the inverse Compton component is sometimes 
complex, because of the contributions of the synchrotron 
self Compton component, especially in the X-ray range, 
and the external component, especially at higher energies. 

As shown in Table 10, the intrinsic luminosities, the 
source dimensions, the bulk Lorentz factors and the view- 
ing angles are quite similar for all sources. The fact that 
RGB J1629+4008 has a steep X-ray (synchrotron) spec- 
trum is explained by our model by the very large value of 
7min, the minimum Lorentz factor of the injected electrons 
which in this case corresponds to 7 p0 ak- The inferred val- 
ues of the magnetic field are not very far from the equipar- 
tition values, with an average ratio {B/B cq ) ~ 1.3. 



7. DISCUSSION 

Our source selection requires the objects to be in the 
HBL region and therefore to have relatively low a ro values. 
Given also their relatively low radio fluxes (and powers; 
see § 7.3), before discussing our results we need to assess 
if our sources are of the same kind as the more power- 
ful and more studied blazars. We do this in the following 
subsection. 

7.1. The Nature of Our Sources 

7.1.1. Are our sources radio-loud? 

Two different definitions of radio-loud quasars have 
been used in the literature (see, e.g., Stocke et al. 1992). 
The first is based on the rest-frame ratio R of radio (typ- 
ically 5 GHz) to optical (typically 4,400 A) flux density, 
logi? = log(/ r /mJy) + 0.4.B — 6.66 (independent of red- 
shift if a r = a Q = 0.5 is assumed). Radio-loud sources 
are taken to have R > 10 which translates, in our nota- 
tion, to a IO > 0.2. The other definition is based on radio 
power and appears to be redshift-dependent, with a di- 
viding line ranging from m 10 24 W Hz -1 for the Palomar 
Green (PG) sample (optically bright and at relatively low- 
redshift) up to 10 26 W Hz -1 for the Large Bright Uuasar 
Survey (LBQS) sample (optically fainter and at higher red- 
shift). The two definitions overlap somewhat (Padovani 
1993). As shown in Fig. 1 and Tab. 1, all our sources 
are firmly in the radio-loud regime as far as the first def- 
inition is concerned. All of them have a ro > 0.26, with 
(a ro ) ~ 0-33 (or (R) ~ 50; typical radio-quiet sources [i.e., 
optically selected quasars] have a ro ~ 0.1 or R ~ 0.5). 
The 5 GHz radio powers of our sources are in the range 
10 24.6 _ 1Q 26 w Hz" 1 , with (L T ) ~ 10 25 1 W Hz" 1 . Given 
their relatively low redshift ((z) ~ 0.2), these values also 
classify our sources as radio-loud. 

7.1.2. Are our sources blazars? 

As described in the Introduction, blazars are character- 
ized by a variety of properties, which have also led to a 
proliferation of names (HPQ, OVV, CDQ, FSRQ). The 
flat-spectrum radio quasar definition is the easiest one to 
apply, and it appears that the majority of powerful FSRQ 
indeed show rapid variability, high polarization, and radio 
structures dominated by compact radio cores, and vice 
versa (Urry & Padovani 1995 and references therein). 

As our sources are less powerful than the FSRQ com- 
monly studied (see § 7.2.4), it is important to check, as 
far as possible, that we are dealing with the same type of 
sources. We do this in the following. 

Our objects have, by selection, a flat radio spectrum. 
In fact, (a r ) = —0.37, that is their typical radio spectrum 
is not only flat but inverted. Three of our sources are 
included in the NVSS and therefore low-resolution (45") 
radio maps at 1.4 GHz are available. RGB J1629+4008 
and RGB J1722+2436 are compact and unresolved. RGB 
J 1629+4008 is also included in the FIRST (Becker et al. 
1995) radio survey, which has a higher resolution (5") than 
the NVSS. This source is still compact and unresolved in 
the FIRST map, with a 1.4 GHz radio flux of 11.9 ± 0.2 
mJy, as compared to 9.0 + 0.5 mJy in the NVSS. The fact 
that the higher-resolution FIRST flux is similar to (actu- 
ally larger than) the NVSS flux suggests that the source is 
core-dominated (modulo variability effects). S5 2116+81 
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has a complex NVSS map, with four radio sources within 
3' from the source position. However, dedicated radio ob- 
servations (Taylor et al. 1996) show that the VLA and 
VLBA (Very Long Baseline Array) 5 GHz core fluxes of the 
stronger source are basically the same, again indicating a 
very compact core. Finally, our own 4.8 and 8.6 GHz Aus- 
tralia Telescope Compact Array (ATCA) maps of WGA 
J0546.6— 6415 (2" and 1" resolution, respectively) indi- 
cate a compact, unresolved source, with a core-dominance 
parameter at 4.8 GHz < 35 (see Landt et al. 2001 for 
details of the ATCA observations). The fact that the 
higher-resolution ATCA flux at 5 GHz is similar to (actu- 
ally larger than) the lower-resolution Parkes-MIT-NRAO 
(PMN) (Wright et al. 1994) flux (339 mJy vs. 287 mJy) 
again suggests that the source is core-dominated (modulo 
variability effects). 

Optical polarization data are available only for two of 
our sources, RGB J1722+2436 and S5 2116+81. The for- 
mer has P opt = 0.70 ± 0.48% (Stockman, Moore, & Angel 
1984) while the latter has P opt = 0.94 ±0.21% (Marcha et 
al. 1996). We note that although these values are below the 
"canonical" division between high-polarized (HPQ) and 
low-polarized (LPQ) quasars (and are indeed not incon- 
sistent with possible scattering by the interstellar medium 
of our own galaxy), many flat-spectrum LPQ are also 
thought to be seen with their jets at small angles w.r.t. 
the line of sight (e.g., Ghiscllini ct al. 1993), and therefore 
classify as "blazars" according to our definition. Further- 
more, if indeed these objects are the equivalent to HBL, 
as discussed in § 7.2, we would expect their duty-cycle 
for high polarization to be lower than their low-frequency 
peaked relatives, as was found for HBL and LBL by Jan- 
nuzi, Smith, & Elston (1994). Finally, the relative promi- 
nent thermal emission, especially in S5 2116+81 (see Fig. 
7), would also explain a low optical polarization. 

As regards variability, at least three of our sources ap- 
pear to vary in the X-ray band, based on our BeppoSAX 
and ROSAT comparison, with RGB J1629+4008 decreas- 
ing in flux by a factor of 4 in 7 hours in our MECS and 
LECS observations. 

7.2. Have we found high-frequency peaked FSRQ? 

The purpose of this paper was to study the X-ray prop- 
erties and, more generally, the SED of a class of flat- 
spectrum radio quasars with relatively large X-ray-to- 
radio flux ratios and effective spectral indices typical of 
HBL. We now address the question if these sources are 
really the strong, broad-lined counterparts of HBL. 

7.2.1. Spectral Energy Distributions 

The SED of RGB J1629+4008 is HBL-like, with the 
steep X-ray spectrum clearly attributable to synchrotron 
emission. For WGA J0546.6-6415 and RGB J1722+2436 
the situation is less clear. Our fits and the SEDs are not 
inconsistent with the fact that the BeppoSAX band might 
be close to the "valley" between synchrotron and inverse 
Compton emission, which would classify them as akin to 
"intermediate" BL Lacs. Finally, S5 2116+81 appears to 
be LBL-like. 

7.2.2. X-ray spectral index and the synchrotron peak 

frequency 



A good indicator of the dominant emission process in 
the X-ray band is the X-ray spectral index - synchrotron 
peak frequency plot (Padovani & Giommi 1996; Lamer, 
Brunner, & Staubert 1996). Padovani & Giommi (1996) 
found a strong anti-correlation between the ROSAT a x 
and fpeak for HBL (i.e., the higher the peak frequency, 
the flatter the spectrum), while basically no correlation 
was found for LBL. This was interpreted as due to the 
tail of the synchrotron component becoming increasingly 
dominant in the ROSAT band as ^ poa k moves closer to the 
X-ray band (see Fig. 7 of Padovani & Giommi 1996). 

The BeppoSAX version of this dependence was studied 
by Padovani et al. (2001; see their Fig. 7). They confirmed 
the ROSAT findings, namely a strong anti-correlation be- 
tween a x and f pca k for HBL and no correlation for LBL, 
with an initial increase in a K going from LBL to HBL (see 
Padovani et al. 2001 for details). 



> 




15 16 17 

Fig. 8. — The BeppoSAX spectral index (0.1 - 10 keV) versus the 
logarithm of the synchrotron peak frequency for our sources (filled 
circles), HBL (crosses), and LBL (open circles). Data for LBL come 
from Padovani et al. (2001), while those for HBL come from Wolter 
et al. (1998; updated in Bcckmann et al. 2002), Beckmann et al. 
(2002), and Padovani et al. (2001). Three of our sources, namely 
those with values of o x < 1, have quite uncertain Vp 0a k values but 
likely in the range 10 14 < i-'peak < 10 16 Hz (see text for details). 



Fig. 8 shows an updated version of the a x — t'peak 
plot for BL Lacs from Padovani ct al. (2001), including 
the HBL studied by Beckmann et al. (2002). We have 
also included in this plot our sources, represented by filled 
points. We notice that RGB J1629+4008 falls right in the 
region occupied by HBL, again confirming our interpreta- 
tion that its BeppoSAX band is dominated by synchrotron 
emission. This is further supported by the detection of 
rapid variability, which has been seen in the synchrotron 
emission of several BL Lacs observed with BeppoSAX but 
never in the inverse Compton component, which varies on 
much longer time scales (e.g, Ravasio et al. (2002) and 
references therein). The other three sources have values 
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which, although quite uncertain, are probably in the range 
10 14 < fp ak < 10 16 Hz, which overlap with the values typ- 
ical of "intermediate" BL Lacs, but display relatively flat 
a x . We notice that the scatter in the diagram is larger 
in the ^ pca k range « 10 14 — 10 15 ' 5 Hz, which is also where 
the HBL/LBL division becomes blurrier. We interpret this 
as due to the fact that in this range of ^ pea k objects can 
be both synchrotron and inverse Compton dominated, de- 
pending on variability state. For comparison, note that 
3C 279, a classical FSRQ, has a BeppoSAX spectral index 
a x = 0.66 ± 0.04 and i/ pcak - 10 13 Hz (Ballo et al. 2002), 
i.e., well into the LBL region. 

7.2.3. The "Blue Bump" Component 

Fig. 7 shows that the thermal component is non 
negligible in most of our sources, with a ratio of black- 
body to synchrotron emission at ~ 5, 000 A ranging from 
- 2 for WGA J0546. 6-6415 and S5 2116+81 to - 0.2 
for RGB J1722+2436. While these ratios are model- 
dependent, especially for the former sources for which the 
optical band is close to the synchrotron cut-off, it is in- 
structive to check the effect of the thermal component on 
the position of our objects in the a ro — a ox plane. By 
subtracting the thermal component and therefore using 
only the synchrotron flux in the optical band to derive 
effective spectral indices, we can see where the objects 
would fall (we assume that the blue bump has no effect 
at 1 keV; but see below). As a rx stays constant while a ro 
and a ox get steeper and flatter respectively, the sources 
move along diagonal lines up and to the left in the di- 
agram. More specifically, the "synchrotron only" values 
are as follows: WGA J0546. 6-6415, a ro - 0.5, a ox ~ 1, 
RGB J1629+4008, a ro - 0.4, a ox ~ 0.85, S5 2116+81, 
a ro ~ 0.4, a ox <~ 1.2 (practically no change is required for 
RGB J1722+2436). It then follows that even by subtract- 
ing the effect of the thermal component our sources would 
still be well within the HBL region in the a ro — a ox plane 
(Fig. 1). Actually, the effect of a strong blue bump, by 
moving sources to the right and down in the diagram, is 
that of "pushing" objects off the HBL region and into the 
radio-quiet zone. Even with some freedom in the mod- 
elling of the synchrotron component, the location of our 
sources in the HBL region results quite robust, and can- 
not be ascribed to the accretion disk. That relatively high 
values of Vpcak are not due to the accretion disk is also 
indicated by the low, "HBL-like" values of a rx , a param- 
eter well correlated with f pea k (e.g., Padovani & Giommi 
1996) and not affected by the optical contamination. Note 
that the v pea ,k values given in Tab. 10 and shown in Fig. 8 
are those relative only to the synchrotron component and 
therefore do not include the thermal component. 

In Seyfert galaxies and radio-quiet quasars the X-ray 
emission is thought to be produced by a hot thermal 
corona sandwiching a relatively cold accretion disk. One 
could then ask if the same component can contribute also 
in our objects, two of which indeed have a relatively large 
disk component in the optical-UV band. While it is cer- 
tainly possible that some fraction of the X-ray emission 
(especially at low energies) is produced by a hot corona, we 
believe that the non-thermal components in these objects 
are dominant, due to the X-ray variability indicated by the 
comparison between BeppoSAX and ROSAT data, par- 



ticularly the rapid LECS/MECS variability seen in RGB 
J 1629+4008, and to the lack of any indication of iron 
line features and hardening of the spectrum due to the 
so-called Compton hump (which are instead common in 
Seyfert galaxies). A possible way to decrease the impor- 
tance of the iron line emission and Compton hump in the 
spectrum could be to assume a dynamic corona, but in 
such case the corona emission would be characterized by 
harder X-ray spectra (a x < 1; see Malzac et al. 2001), and 
thus could not account for soft excesses or steep spectra. 

7.2.4. Source Powers 

As discussed above, the 5 GHz radio powers of our 
sources are in the range 10 24 6 - 10 26 W Hz -1 , with 
(L T ) ~ 10 25 - 1 W Hz" 1 . Their 1 keV X-ray powers are 
in the range 10 19 ' 4 - 10 20 6 W Hz" 1 , with (L x ) - 10 20 1 W 
Hz -1 . These values, although not very different from the 
mean values of the EMSS BL Lacs ((L r ) ~ 10 248 W Hz" 1 
and (L x ) - 10 20 ' 3 W Hz" 1 ), are still in the FSRQ range, 
albeit on the low side (see, for example, Fig. 5 of Perlman 
et al. 1998 and Padovani et al., in preparation). Their 
radio powers put our sources near the low-luminosity end 
of the DXRBS FSRQ radio luminosity function (Padovani 
2001). 

7.3. Astrophysical Interpretation 

We have found three examples of strong, broad-lined 
blazars with an SED typical of HBL (one source) and in- 
termediate BL Lacs (two sources). While the first case is 
quite strong, the other two are more tentative. These are 
the first such sources, as most broad-lined blazars known 
so far have SED typical of LBL, i.e., with synchrotron 
peaks at IR/optical energies. 

Our FSRQ, however, have relatively low powers, more 
typical of BL Lacs than quasars. This ties in with the sug- 
gestion by Ghisellini et al. (1998) and Ghisellini, Celotti, 

6 Costamante (2002) that blazars form a sequence, con- 
trolled mainly by their bolomctric luminosity which in turn 
controls the amount of cooling suffered by the emitting 
electrons. In powerful sources, where cooling is severe 
(high radiative (U r ) plus magnetic field (Ub) energy den- 
sities), electrons of all energies cool rapidly (i.e., in a time 
shorter than the injection time), making the particle dis- 
tribution to have a break at 7 = 7 m i n , which in this case 
becomes 7 pea k- These sources are typically LBL-like, and 
are those located in the lower part of the 7 pca k — U T + Ub 
plot (see Fig. 9). In less powerful sources, instead, only the 
highest energy electrons can cool in the injection time, and 
consequently the particle distribution will have a break at 

7 = 7c, for which i coo i = i inj . In this case 7 pca k = 7c (and 
oc (U T + Ub)^ 1 ', see details in Ghisellini, Celotti, & Costa- 
mante (2002)). Fig. 9 shows that the FSRQ studied in 
this paper, which have relatively high 7 pca k (and therefore 
Vpeak) values, appear to fit within this scenario, as they 
follow the suggested sequence. (Note that, given the scat- 
ter in the plot, this conclusion stands even in the case of a 
decrease of <~ one order of magnitude in the 7 pca k values 
for the three sources with relatively uncertain estimates of 
this parameter.) Moreover, in agreement with their rela- 
tively low powers, they are positioned at the low-end of 
the FSRQ region and actually within the BL Lac region. 
In other words, we still have not found a case of high-f pea k 
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(and therefore high-7 pea k) and high-power blazar, which 
would invalidate this scenario. 
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Fig. 9. — 7 pca ki the Lorentz factor of the electrons emitting most 
of the radiation, vs. the total energy density in the emitting re- 
gion, i.e., the sum of the radiative (U r ) and magnetic field (Ub) 
energy densities, as measured in the comoving frame. Open circles 
represent "extreme" HBL, filled circles represent TeV sources, open 
squares represent BL Lacs, open triangles represent FSRQ, filled 
triangles are FSRQ at z > 4, while stars indicate our sources (HF- 
SRQ). All data apart from the latter are from Ghisellini, Celotti, & 
Costamantc (2002), updated in part from Ghisellini et al. (1998). 

Another explanation for the relatively low power of our 
sources is related to a selection effect. Fig. 7 shows that 
our only "true" HBL-like FSRQ, RGB J1629+4008, has 
a synchrotron peak at ~ 2 x 10 16 Hz, well above the peak 
of the thermal emission. Its non-thermal power is also low 
enough (as implied by its radio power Lsghz ~ 6 x 10 24 
W Hz -1 ) that the thermal component, and its associ- 
ated broad lines, is non-negligible in the optical/UV band. 
Given that the equivalent width of its strongest emission 
line is ~ 80 A, the same object with a non-thermal power 
a factor of ~ 20 higher, corresponding to Lsqhz ~ 10 26 W 
Hz -1 , still rather modest for an FSRQ, would be classi- 
fied as a BL Lac object, since its equivalent width would 
go down to ~ 4 A. An even stronger non-thermal compo- 
nent would completely wash out the emission lines, mak- 
ing the redshift determination of this object impossible, 
as there are no galaxy absorption features in its spectrum. 
It then follows that a broad-lined source with the same 
SED as RGB J1629+4008 but a much larger power would 
not even be classified as an FSRQ but as an HBL without 
redshift. This could also explain why high-power, high- 
^poak BL Lacs have not been yet found. A key ingredient 
of this argument is the frequency of the synchrotron peak. 
Assuming a typical ^ pca k value for a "classical" FSRQ of 
~ 10 13 ' 5 Hz, which implies a much weaker optical syn- 
chrotron component, an object with a radio power ~ 100 



times larger than RGB J1629+4008 will have the same 
optical (non-thermal) power and therefore will still be rec- 
ognized as a quasar (assuming similar thermal powers). It 
follows that high-power, low-f pea k FSRQ can be identified 
as quasars up to very high luminosities. Conversely, high- 
t'peak blazars would be classified as HBL without redshift 
already at moderate radio luminosities. One problem with 
this explanation for the lack of high— v pe ak — high-power 
blazars is that there should be intermediate cases, i.e., 
HBL with broad (but weak) lines. This does not seem ap- 
parent, for example, from the optical spectra of the EMSS 
BL Lacs (e.g., Rector et al. 2000). 

As the selection of our sources was done 3 — 4 years 
ago, one could ask if our sources are still really extreme 
in terms of their X-ray-to-radio flux ratios. Out of the 
199 DXRBS FSRQ known to date, WGA J0546.6-6415 
is still the second most extreme. The FSRQ RGB sam- 
ple, given its higher X-ray flux limit and slightly lower 
radio flux limit, is even better suited to find sources with 
high / x //r ratios (lower Padovani et al., in prepa- 

ration). Sorting the known RGB FSRQ in order of as- 
cending a rx value, RGB J1629+4008 is at the top of 
the list while RGB J1722+2436 is number four (RGB 
J1413+4339, which is number two, is a borderline quasar; 
see Laurent-Muehleisen et al. 1998). Finally, S5 2116+81 
is number 15 out of the ~ 400 FSRQ in the multifrequency 
catalog of Padovani, Giommi, & Fiore (1997) with X-ray 
data sorted in order of ascending a rx value. It was selected 
because it is the X-ray brightest. 

We find it significant that, having selected some of the 
few known FSRQ with effective spectral indices typical 
of HBL and therefore relatively large X-ray-to-radio flux 
ratios, we found only one case where the X-ray band 
is clearly dominated by synchrotron emission, with a 
^peak ~ 2 x 10 16 Hz, while the other three sources have 
(z-'peak) ^ 10 15 Hz. As shown, for example, in Fig. 8, HBL 
can reach much higher values, up to 10 18 — 10 19 Hz, albeit 
in cases of extreme variability (e.g., MKN 501; Pian et al. 
1998). A more meaningful comparison is with the HBL in 
the multifrequency catalog of Padovani, Giommi, & Fiore 
(1997) having radio flux > 20 mJy, the minimum value 
for our sources. In that case the u pca ^ distribution has a 
typical value ~4x 10 16 Hz and reaches ~ 2 x 10 17 Hz. 
A discussion of the ^ pca k distribution for FSRQ and BL 
Lacs belonging to the same sample will be presented by 
Padovani et al. (in preparation). 

Sambruna, Chou, & Urry (2000) reported on ASCA ob- 
servations of four FSRQ characterized by steep ROSAT 
spectra (a x ~ 1.3). The sources were all found to have 
flat hard X-ray spectra, with a x ~ 0.8. Sambruna et 
al. discuss their results in terms of relatively high syn- 
chrotron peaks and thermal emission extending into the 
X-ray band. We stress, however, that their sources sam- 
ple a region of parameter space widely different from ours. 
Their effective spectral indices, in fact, place their FSRQ 
firmly in the LBL region, unlike ours, so that their four 
objects should not have been expected to show high v pea k 
values and steep ASCA spectra. Georganopoulos (2000) 
used the fact that these four FSRQ happen to have rela- 
tively low core-dominance parameters (« 1) compared to 
two "intermediate" BL Lacs to suggest that these sources 
could be seen at larger angles than more typical FSRQ. 
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We stress that, as discussed in § 7.1.2, our sources, and 
especially RGB J1629+4008, appear to be quite compact 
so that his interpretation might not apply to our objects. 

8. SUMMARY AND CONCLUSIONS 

We have presented new BeppoSAX observations of four 
flat-spectrum radio quasars selected to have effective spec- 
tral indices (a ro , a ox ) typical of high-energy peaked BL 
Lacs. The main purpose of the paper was to see if these 
sources are indeed the broad- lined counterparts of BL Lacs 
with the synchrotron peak at UV/X-ray energies (HBL). 
Our objects are quite extreme in terms of their X-ray- 
to-radio flux ratios (a factor ~ 70 higher than "classical" 
FSRQ), qualify as radio-loud sources both in terms of their 
radio-to-optical flux ratio and power, and are variable and 
radio-compact. Our main results can be summarized as 
follows: 

1. We have discovered the first FSRQ (RGB 
J1629+4008) whose X-ray emission is dominated 
by synchrotron radiation, as clearly shown by our 
BeppoSAX observations (a x <~ 1.5 in the 0.1 — 10 
keV band), in which the source was also rapidly 
variable. This object is therefore the first confirmed 
member of a newly established class of HBL-like 
FSRQ. This result is fully consistent with earlier 
ROSAT observations which detected a steep soft X- 
ray spectrum (a x ~ 2 in the 0.1 — 2.4 keV band). 
The combination of the X-ray data with archival 
radio and optical data gives a spectral energy distri- 
bution typical of HBL sources, with a synchrotron 
fpeak ~2x 10 16 Hz. The derived values follow very 
well the a x — f pca k correlation seen for HBL sources 
(Fig. 8). 

2. The other three sources show a relatively flat (a x ~ 
0.75) BeppoSAX spectrum. However, we have 
found for two of them (WGA J0546-6415 and RGB 
J1722+2436) some indication of steepening at low- 
encrgics. This is based on BeppoSAX and ROSAT 
data and on their (sparsely sampled) SEDs. We 
interpret this as the tail of synchrotron emission, 
which is not strongly constrained but would peak 
around « 10 15 Hz, as typical of "intermediate" BL 
Lacs for which the synchrotron and inverse Comp- 
ton components meet within the BeppoSAX band. 
The last source (S5 2116+81) turned out to be a 
more typical FSRQ, with the BeppoSAX band fully 
dominated by inverse Compton emission and a syn- 
chrotron ^pcak < 10 15 Hz. 

3. By fitting a synchrotron inverse Compton model 
which includes the contribution of an accretion disk, 
whose power we estimate from the broad-line region 
luminosity, to the spectral energy distributions, we 
have derived physical parameters (e.g., intrinsic 
power, magnetic field, etc.) for our sources. The 
thermal component was found to be non-negligible 



in three out of four objects, with ratios of thermal 
to synchrotron emission in the range 0.6 — 2 at 5,000 

A. 

4. Although the original selection was based mostly on 
the effective broad-band (radio through X-ray) col- 
ors, which are in first approximation independent 
of luminosity, all four candidates are at relatively 
low powers, more typical of BL Lacs than of FSRQ 
and in any case close to the low-luminosity end 
of the FSRQ radio luminosity function (Padovani 
2001). Wc interpret this as due to two, possibly 
concurrent, causes: an inverse dependence of 7 pca k 
(and therefore ^ pca k oc 7p 0ak <5-B), the Lorentz fac- 
tor of the electrons emitting most of the radiation, 
on bolometric power, due to the more sever cooling 
at work in more powerful sources; and a selection 
effect, namely the fact that a high-power, HBL 
like quasar will have such a strong optical/UV non- 
thermal component that its emission lines will be 
swamped and the object will be classified as a fea- 
tureless BL Lac. 

A better understanding of these rare sources and the 
confirmation of the non-existence of high-^ pca k — high- 
power blazars will require the selection and study of a 
larger sample of objects. 
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